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Low energy string theory predicts the existence of an axion field which can lead to cosmic 
birefringence. We solve the electromagnetic wave equations in the presence of such an axion and a 
dilaton field in order to determine their effect on the polarization of light. We find that the presence 
of dilaton field leads to a nontrivial modification of the final result. We comment on the possibility 
of discovering such an effect by observations of radio wave polarizations from distant radio galaxies 
and quasars. We have also determined the limits on the string theory parameters that are imposed 

. . . by the current radio polarization data. 
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• I. INTRODUCTION 

(D ; 

' There have been several studies which explore the possible existence of cosmic birefringence ||] . In recent work Kar 

, et al Q and Majumdar and Sengupta ^ has pointed out that the Kalb-Ramond field, which arises in supergravity 

' theories, makes the space time birefringent. The authors solve the electromagnetic wave equations in the presence of 

l such a field and predict that this field will lead to a redshift dependent rotation of polarization whose details depend 



> 



on whether the universe is radiation or matter dominated. Earlier in Ref. [Q Carroll et al analyzed the polarizations 
of radio waves from distance galaxies and quasars as a function of redshifts and found no statistically significant effect. 
[ — r The observable of interest in that study was the angle f3 = x — ip, where -0 is the orientation angle of the axis of the 
' radio galaxy and x the observed polarization angle after the effect of Faraday rotation is taken out of the data. This 
can be done by making a straight line fit to the polarization angle ^(A^) which can be expressed as 



o 
o 



0{X^) = RMX^+X (1) 



' where A is the wavelength, RM and x a-re the slope and the intercept of the fit. RM stands for rotation measure and 
depends on the plasma density and the parallel component of the background magnetic field along the direction of 
propagation of the wave. Later Ref. Q claimed the existence of a large scale anisotropy in the same data. This claim 
I was, however, questioned by many authors j^^. 
f~i . The possibility of local anisotropy, independent of redshift, was explored in Ref. pland car her in (1O0. These 
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authors found a statistically significant effect. The original claim made by Birch 1 101] in 1982 was dismissed by 
Bietenholz and Kronberg |l2[ in 1984 who compiled a larger data set which they claimed did not show any effect. 
However Jain and Ralston argued that Ref. did not pay careful attention to the parity of the correlation 
M I ansatz [ p^ . When the transformation properties of statistical distribution are taken into account the data does show 
- - I a significant effect which is further enhanced by choosing an appropriate statistic as well as by putting some cuts on 
the data. 

Inspite of the progress in our understanding in string theory in the last two decades, we are yet to have a string 
theory prediction that can be tested in the experiments. One of the main reason for such an unsatisfactory situation 
is that the energy scale at which stringy effects become significant is beyond the reach of present day experiments. 
It is thus believed that the most likely area for a confrontation between string theory and experiments is through 
extracting astrophysical or cosmological consequences of string theory which might be measurable today. With this 
kind of motivation in mind, we study in this paper the prediction of string theory for cosmic birefringence. However, 
as our knowledge of string theory is mostly perturbative in nature, we face several problems in carrying out such a 
study. These are mainly due to the presence of various massless moduli that are present in the low energy limit of 
string theory. We will point those out as we proceed. 
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Four dimensional heterotic string theory, in its low energy limit, contains moduli coming from compactification, 
gauge fields, dilaton and an axion. In this theory, coupling among various fields are unique and they are all controlled 
by the expectation value of dilaton. Given such a framework, we analyse the propagation of electromagnetic waves in 
curved metric background. Due to the presence of dilaton, as we will argue, medium turns into a conducting medium 
where the time variation of dilaton plays the role of effective conductivity. On the other hand, the axion rotates the 
plane of polarization of the electromagnetic wave causing non-trivial birefringence. However, notice that the dilaton 
(and other moduli) is a massless scalar in string theory. An important issue to resolve in string cosmology is to explain 
the absence of massless dilaton at the present time. One expects that a potential for dilaton will be generated in string 
theory through some non-perturbative effects and the dilaton will sit at the minimum of the potential picking up a 
mass. The mechanism of how precisely this should happen is still unknown in string theory. This, in turn, creates 
the most severe problem in comparing observational data with string theory. In the rest of the paper, we continue 
to treat dilaton as a massless scalar and find string prediction for cosmic birefringence. We hope that some of our 
results will be helpful when dilaton decoupling in string theory is properly understood. In general, we find that the 
relevant equations are hard to solve exactly. However, it is well-known that heterotic string theory in four dimensions 
is self dual. The duality group is parametrized by SL{2,R) matrices. Using this property, in the next section, we 
generate some non-trivial solutions for various fields starting from some simple seed configuration. This, in turn, 
allows us to calculate birefringence of electromagnetic waves, propagating in radiation or matter dominated universe, 
in a fairly straightforward manner. However, since the solution generating technique does not saturate all possible 
allowed configurations, in section HI, we solve the electromagnetic equations within the WKB approximation. We 
then discuss our results in section IV. 



II. BIREFRINGENCE FROM HETEROTIC STRING 

In the low energy limit, the heterotic string compactified on six-torus is represented by the following action [l4| : 
S = J d^xV^e-'^ [Rg + G^'dp^d.cf) - j^Hp.pH^'^P 

~ G^'PG^'^Fl^J (LML)abF^''\x + ^G^^^Tiid^MLd^ML)] (2) 
where the gauge fields and the anti-symmetric three rank tensor field are defined as 

Hpvp = {dpB^p + 2A^i^'^LabFl^J) + cyclic permutations of /i, v, p. (3) 

The superscript a on runs from 1 to 28 denoting the 28 gauge fields that appear in the four dimensional heterotic 
string action. Out of these 28 gauge fields, 16 are originally present in ten dimensions. The rest 12 appear due to 
torus compactification. Among other terms in (||), Rq is the scalar curvature associated with the metric Gp^ and 
is the dilaton. The 28 x 28 dimensional matrix M encodes all the scalars that appear due to compactification. Q In 
what follows, the detail of their structure will not be important. The matrix L in (||) is also 28 x 28 dimensional and 
is given by 

/O /g \ 

i = /6 (4) 
V -/i6 / 

where /„ denotes n x n identity matrix. 

In the following, however, we will work with a truncated version of the effective action. We will turn off all the scalars 
that appear due to compactification. This essentially amounts to setting the matrix M to identity. Furthermore, we 
will set all the gauge fields to zero except one. We will take this one to be present in the original ten dimensional 



^ Note that in (g), the first two terms in the action is like Brans-Dicke (BD) theory with the BD parameter lo = ~1. The 
radar-echo delay experiments jlHj in solar system have set the limit on u > 500. So in a sense, the theory as it stands is ruled 
out by this experiment. However, let us note that the string action in general has many moduli field (such as matrix M here). 
If these scalars happen to have flat direction along the dilaton, the effective lj might change from its value —1. 
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heterotic action. We will see that this choice of trunctation makes the dynamics tractable and at the same time keeps 
the essential physics intact. In this approximation, the action, in the Einstein frame, reduces to 

J d^xV^lRg - \crd,ct^dA - l-e-^^H^.pH^-P - e^^F.^Fn (5) 

The equation of motion for the S^i, following from the above action is 

Df,{e^^'f'H''''P) = 0. (6) 
This, in turn, allows us to introduce a pseudoscalar field known as axion, through the relation 

TJM^P „_l^e2'^e^'''"'9^«'. (7) 
With this definition of the rest of the equations of motion are the following 

= ^e^^'d^m.'f + 2e-^F^pFP - Ig^.e-^Fp^F^- + ^d^^d,<P, (8) 

Dp{e-'''F^"' + -^F^"') = Q, (9) 

g20^M£,^^ + 2e^'''d^,(t>d^'^ ~ Fp^.F'"' = 0, (10) 

e^'^D'^Dpie-''') + e^'/'d^^d"^ - e'*9^(/.9''0 - Ff.^F"" = 0. (11) 

Here 

i^^" = -^e'"''"^Fp,, (12) 



which satisfies Bianchi identity 

DpF'''' = 0. (13) 

It is easy to check that the set of equations (|l)-(0) and ( p^ ) are invariant under the following set of SL(2, R) 
transformations [Q 

A^A' = ^^, F^,^F'^^ = {c^ + d)Fp,-ce-*Fp,, g^. ^ g^., (14) 

where A = 4" + ie~'^ and a, b, c, d are real parameters with ad — he = 1. 

Here, in what follows, we will restrict ourselves to the case where dilaton, axion and the metric are only function 
of time. Our startegy will be to solve all the equations of motion in a fixed metric background. First, we take the 
metric as 

ds^ = R'^{ri){-drf + dx^ + dy^ + dz^). (15) 
The gauge field strength in this background can be defined as fl^ 



I —Ere —Ey —E. 

Ex Bz —B, 

Ey —Bz Bx 

\Ez By —Bx 



(16) 



With this definition, equation (^ reduces to 

V-E = 0, a„(e-'^Ei?2)-Vx (e-"^Bi?2) = -2a„*Bi?2 (17) 
Similarly the scalar equations are 

^"^"^ " J?^''^^"* ' -^^Mdv-^ + 4i?'e-2'^E -6 = 0, (18) 

2 p^"^ 

a^aV - j^drjRd^c^ + -j^{d,,^ f - 2R\-1'{E^ - B^) = 0. (19) 
On the otherhand, ( p^ ) reduces to 

V • B = 0, d^{R^B) + V X (i?2E) = 0. (20) 
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A. Analyzing Equations in flat background 



Before we go on to analyse more realistic scenario, it is worthwile to study ( |18| ) - (|2y) in flat metric background for 
which R{ri) — 1. Beside being completely tractable, this will indeed give us an idea of dilatonic and axionic influence 
on the propagation of electromagnetic waves. In this very special context, equations (|l^)-(|2^) reduce to the following 
forms: 



d^cj) - e^'I'idrj-^f + 2e-'^(E2 - B^) = 0, 
a^(e-^a^B) - V2(e-'^B) = 2a„*V X B. 



(21) 
(22) 
(23) 



It is easy to check that when (f> is turned of, the above equations reduce to the case analyzed in |Q] . The very presence 
of dilaton (f) effects the dynamics in a very crucial way. In order to understand the behaviour of the solutions, we first 
set the axion field to zero. Then equation (ESh reduces to 



dp 



drjCfldrj'B 



V^B = 0. 



(24) 



Such an equation for magnetic field is quite well known from the studies of electromagnetic wave propagation in 
conducting medium. As we can see from (^), the time variation of dilaton plays the role of effective conductivity of 
the medium. Such a damping factor, as usual, reduces the amplitude of the magnetic field B. 

The set of equations (p0|)-(p3|) are in general very hard to solve directly. Though it would be desirable to have such 
solutions, we would however use the property ( p^ ) of heterotic string theory to generate new solutions from some 
simple ones. To this end, we notice that if we set the dilaton field to zero (using sift symmetry of dilaton in heterotic 
string theory, we can always bring any constant dilaton field to zero), equations (p0[)-(p3|) possess a set of solutions 0] 



* = hdt = ht, B±{t, z) = b±e 



ikz ^iLJ± t 



(25) 



where h is an integration constant and uj± — v«?±2/ifc. In the above equations, we have defined B(t, z) — b(t)e*'^^ 
with z being the direction of propagation. Furthermore h± are defined as the circular polarization states, namely 
b± = hx^ iby. Now writing B± = &-|-e**^^e*'^± and birefringence as Sip = \{^+ — (fi-), we get 



Sifi = h(t2 — <i), for k >> h. 



(26) 



Here, for flat metric, rj is same as cosmic time t. Thus in the above equation, ti and t2 are the earlier and present 
time. 

As discussed earlier in (|l^, heterotic string theory in four dimension has self-duality symmetry. As a consequence, 
equations of motion for various fields are invariant under SL{2,R) transformation. These transformations have the 
property to generate non-trivial field configurations from the known ones. A generic solution is thus parametrized 
by three independent parameters of SL(2,R) group. Here, we start with the field configuration given in ( p5| ) where 
dilaton field is set to zero. Next, by the use of SL(2, R) symmetry, we generate configurations with nontrivial dilaton 
and various other fields. Strictly following the transformation rules given in (p^, and denoting the new configurations 
by (*',0',B',E'), we, thus, get 



ac*^ + (be + ad)^' 



bd) 



.-0' - 



1 



(c* + dY + c2 



(c* + d)2 + 



(27) 



where ^ is given in (p5|). As expected the solutions are parametrized by three independent SL(2, i?) parameters 
(taken here as) 6, c, d. The other parameter a is determined from the relation ad — be = 1. Furthermore, from the 
transformation property of F^j_i, under SL{2, R), we get 



B'± = [(c* + d) T ^c]b±. 



for k » h 



Thus the birefringence is then given by 

V = h{t2 - h) 



tan 



_i / 2(c*-|-o?)c 
\{c^ + d)^-c^ 



for fc >> h. 



(28) 



(29) 



4 



B. Analyzing the equations in curved space 

Here, we will study (H)-(|3|) for the general metric given in (^5|). We will follow the same strategy as the previous 
subsection in order to find solutions of equations of motion. Keeping this in mind, let us first consider (|^)-(1^) when 
the dilaton is set to zero. In that case, defining F± — B?{'q)B±, we get an equation for magnetic field as follows ||^: 



dr], 



(30) 



with h being a constant. If R{'ri) increases fast enough with rj (which is indeed the case for radiation and matter 
dominated universe), a WKB type analysis will give us a very close to exact result for the above equation. Carrying 
out such a computation we get 



(31) 



where we have S± = ± J 



Like B± as in eqn.(24), one can also rewrite the eqn.(30) in the following way, 

F± = /±e*'=^e*"±'' (32) 



where, uj±_ 



Now we have, 
B± 



-F+ 



and E± = ± 



kB?(r\) dr] 



(33) 



From the above experssions for B± , we see that the birefringence due to the presence of axion (dilaton is set to zero) 
is given by. 



ni 



hdt 



(34) 



where we have converted the conformal time (77) into cosmic time t and ti and ^2 are the earlier time and present 
time respectively. 

As before, starting from ( |30| ) and ( p3| ) and using the freedom of SL(2, R) symmetry of equations of motion, we 
generate here large class of configurations, where axion and dilaton are given by 



ac^'^ + {be + ad)* + (ac + bd) 



1 



(C* + d)2 + c2 ' 



{Cf + df + c2 

where ^ is given in (|30|). On the other hand, from the transformation property of Ffj_i, under SL{2, R), we get 

h 



B' 



(c* + d) T «c(l ± 



B±. 



Consequently, the expression of birefringence is given by 



6ip' = — tan ^ 
^ 2 



2(0^- + d)c 



(Cf +d)2-c2(l-^) 



')1 



''^^ hdr] 



(35) 



(36) 



(37) 



After having 5lp' for general metric, let us concentrate on two particularly interesting cases; first radiation dominated 
case and then matter dominated case. 

Bl. Radiation dominated Universe: For radiation dominated universe, R is given by R{t) — RoRt^^^, where Rqr 
is a dimensionful constant. Instead of writing various fields explicitly in this particular case, we give the expression 
of birefringence that follows from (|36| ) . It is given by 0: 



^The relation that we used between R{t) and z can be described as follows, 
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1 

- — tan 
2 



2c(c* + d) 



(39) 



B2. Matter dominated Universe: For matter dominated phase, R is given by R{t) = Rmit^''^- In the same way 
as before, one can show that the birefringence is now given by: 

dip = — tan 



where, in both cases, the second term corresponds to the birefringence corresponding to zero dilaton field 0. 
However, in both the cases, we see that the expression of <5(p' depends on redshift ^ in a crucial way. Furthermore, 
though very weak, has a dependence on the frequency k of the electromagnetic wave. 



2c(c* + d) 



l-(l + z) 



3/2 



(40) 



III. APPROXIMATE SOLUTIONS INCLUDING THE DILATON FIELD 



So far we have concentrated on solving the equations by first setting the dilaton field equal to zero and then 
using the SL(2, R) transformation to generate solutions with nonzero dilaton field. However this procedure does not 
guarantee that all possible solutions to the equations are generated. We next solve the field equations directly within 
the WKB approximation. We first write the electromagnetic field equation in terms of the variable S±. Since S± is 
assumed to be a slowly varying function we drop its second derivatives and terms proportional to [iSj-]^. With these 
approximations the field equations reduce to 



'dt 

Si 



^R.'dt'dt 

-4>I2 



R^' 



-20 



(41) 
(42) 

(43) 
(44) 



where and 5*4 are the real and imaginary parts of the function S± . We also drop the term proportional to in 
the dilaton field equation (|4^) since it is higher order in h and hence expected to give negligible contribution. We 
then find that Eq. 4^ can be solved to give 

(t)^ B + A/t 



where A and B are constants. 



Furthermore, are given by. 



S^ 



he 



ARIj^j 



1 — exp 



(45) 



The rotation in polarization is given by 5(f) = {S+ — S-)/2. The functional form of the rotation angle is plotted in Fig. 
1 as a function of redshift for some arbitrary choice of the parameters ( = —e^^/*''^he^^/AR^j^,j and ^ = A/t2- Hence 
we find that the dilaton field leads to a nontrivial change in the redshift dependence of the rotation in polarization. 
Whereas in the absence of dilaton we predict a steady increase of the polarization rotation angle with redshift, in the 
present case the rotation turns on only at large redshifts. 



#4 = 1+,. (38) 

R{ti) (ti)" ^ ' 



where for Radiation dominated phase a = \ and for matter dominated phase a = | 
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FIG. 1. The redshift dependence of the polarization angle for some representative values of the parameters ^ and ^ 

We next consider the polarization data from distant quasars and radio galaxies in order to determine if there exists 
a signal of the type predicted by our analysis. We considered all the sources that were used in the statistical analysis 
in Ref. |^ for which the information about redshift was available. The NASA extragalactic data base (NED) was used 
to update the information about redshifts. We compiled a total of 231 sources in this manner^. We also examined 
the data after removing all sources which lie within ±30° of the galactic plane. This cut is useful since it removes 
many sources with very large rotation measures and hence reduces the possibility of bias in data. The bias in rotation 
measure can arise due to the mr ambiquity in the measurement of linear polarization. After imposing this cut on 
the data we are left with a total of 160 sources. We use the maximum likelihood analysis in order to determine the 
presence of correlation in data. For the null hypothesis we use the von Mises (vM) distribution which serves as a 
prototype for the statistical fluctuations for circular data. It is given by |[l7[|, 



/vm(/3) 



exp [kcos2(/3 — /?)] 
7r/o(K) 



(46) 



where k is a parameter which measures the concentration of the population, (3 is the angle between the axis of the 
radio galaxy and the observed polarization after taking out the effect of Faraday rotation and /3 is the mean angle. 
The factor 2 has been inserted since the polarization angle is ambiguous by mr and not 2n7r JTst . The distribution 
peaks very close to /? = 7r/2 and therefore it is reasonable to eliminate this parameter by setting it exactly equal to 
tt/2 We are therefore left with only one free parameter and the fit to the entire data set leads to k = 0.67 and 
after the galactic cut we find k = 0.71. 



^The entire data set is available on the website home.iitk.ac.in/~pkjc 
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FIG. 2. The distribution of the offset angle /3 between the intrinsic polarization x and the orientation axis of the radio galaxy 
^ for (a) full data containing 231 points and (b) after putting the cut on redshift, z > 0.5. 

We next study the redshift, z, dependence of distribution. It turns out that the distribution has a strong dependence 
on the redshift which can be modelled by the following generalization of the vM distribution, 



z) = N{k, Ki)exp [-{k + kiz) cos(2/?)] 



(47) 



where N{k, Ki) is a normalization factor. Maximizing the likelihood in this case gives k = 0.29 and ki — 0.90. The 
maximum likelihood increases by 7 units, i.e. Ljz — Lj = T = 7, where L/ and L77 are the maximum likelihoods 
for the distributions given in equations and ( p7| ) respectively. After making the galactic cut the corresponding 
parameters are k = 0.32, ki = 0.97 and T = 5.1. The statistic 2T is distributed as Xi- The increase in the maximum 
likelihood is quite large and implies that the width of the distribution has significant dependence on the redshift. This 
is also apparent from a plot of the distribution of the angle 2/3 shown in Fig. ^ for the complete data set and after 
imposing the cut on redshift, z > 0.5. We find that after imposing the cut the distribution is more strongly peaked 
at /3 = 7r/2. 

In order to search for a correlation of the type implied by Eq. we propose the following distribution 



G{f3, z) = iV(K, Ki, C, Oexp [-(k + Kiz) cos(2/3) + C, sin(2/3)] 
where iV(K, ki, C, ^) is a normalization factor and C^, ^) is given by 



1 — exp 



3/2 



(48) 



(49) 



This form of the distribution will lead to the relationship ( |45| ) for the polarization angle < /3(z) >, where <> denotes 
the mean value, as long as deviation of < /3 > from 7r/2 is small. 
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The best fit parameters for the correlated ansatz are k = 0.27, ki = 0.98, ( = 8.0 x 10 , ^ = 7.1 with the statistic 
2T = 2{Liii - Lii) = 5.34 for the full data set and k = 0.28, m = 1.1, C = 8.8 x 10"^, $ = 4.1 with the statistic 
2T = 6.9 for the case of galactic cut. Here Lju is the maximum likelihood for the distribution given in equation (|4^). 
The corresponding percentage values (p-values) for these two cases, with and without the galactic cut, are 7% and 
3.5% respectively. These p-values represent the probabilities that the correlation seen in data can be obtained from a 
statistical fluctuation. After making the galactic cut the results correspond to a 2(T effect. Hence these results show 
marginal statistical significance. Although they show a significant correlation it is clear from the best fit parameters 
that the solution prefers the extreme values of the parameter C, very close to zero. This implies that the deviation of 
(3 from 7r/2 is very close to zero for a wide range of redshifts and deviates from zero only for relatively large redshifts. 
This is clearly seen in Fig. 3 where we plot the data and our best fit model for the case of galactic cut. On the y-axis 
we plot the sin 2/3 in order to deal with quantities invariant under angular coordinate transformations. The data, 
sin 2/3, has been averaged over the redshift interval of 0.25 in order to reduce the noise. The number of data points in 
each bin, starting from the smallest redshift, are 54,19,15,13,10,9,9,11,7,4,3,2,2,1,1 respectively. 




0.5 1 1.5 2 

z (redshift) 



FIG. 3. The fit to the polarization offset angle /3 for distant radio galaxies and quasars. The redshift is plotted on the x-axis 
and the mean (over a redshift interval of 0.25) of sin(2/3) is plotted on the y-axis. The region within ±30° of the galactic plane 
was eliminated from the data. 

If we eliminate the two largest redshift points then we do not see any significant effect. Strictly speaking we were 
unable to find the maxima of the likelihood. The fit was driven to ridicuously large values of ^ > 100 and we did 
not search beyond these extreme values. The statistic 2T in this case was roughly 1.6 which has p-values of 45% and 
hence shows no significant correlation. In Fig. 4 we show the limits of the parameters C and ^ that are imposed by 
the polarization data. The limits are obtained by excluding all sources which lie within ±30° from the galactic plane. 
However the results remain unchanged even if this cut is not imposed. 
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FIG. 4. Limits on the values of the parameters and ^ from the polarization data from distant quasars and radio galaxies 
after eliminating all sources which lie within ±30° from the galactic plane. The results do not change appreciably if this cut is 
not imposed. The shaded region represents the allowed range of these parameters. The two largest redshift points were deleted 
for this analysis. 

IV. CONCLUSIONS 

In this paper, we have studied cosmic birefringence in the framework of perturbative four-dimensional heterotic 
string theory. Due to coupling between dilaton and axion in this theory, we have shown that the dilaton field changes 
the birefringence produced by the axion field in a nontrivial way. Wc obtained an approximate expression for the 
rotation in polarization predicted by heterotic string theory within the WKB approximation. Wc analysed the data 
for distant quasars and radio galaxies to determine if the effect is present in the data. We find a marginally significant 
effect. However the correlation is caused primarily by two largest redshift points and the effect is lost after these 
points are eliminated. Hence we find that the data does not show significant correlation but the large redshift points 
are suggestive that a signal might exist. This should be thoroughly investigated by accummulating more data at the 
largest redshifts. We have also determined the limits on the string theory parameters that are imposed by the current 
polarization data from distant quasars and radio galaxies. 

Wc, furthermore, like to point out that all these calculations were done taking dilaton as a masslcss scalar. As 
explained in the introduction, one of the biggest problem in string theory is to find a mechanism to remove this 
massless dilaton from present day physics. One expects that a potential for dilaton will be generated in string theory 
through some non-perturbative effects and the dilaton will sit at the minimum of the potential picking up a mass. 
As the mecanism of how precisely this should happen is still not clear in string theory, we have not dealt with such a 
scenario in this paper. 
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